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ABSTRACT 

We  present  a  new  approach  toward  robust 
signal  detection  which  is  based  on  techni  - 
ques  rooted  in  differential  geometry. 

These  methods,  as  opposed  to  the  commonly 
employed  classical  saddlepoint  criteria, 
readily  admit  the  quantitive  measure  of  the 
degree  of  robustness  over  very  general 
classes  of  admissable  noise  distributions. 
Our  approach  thus  is  seen  to  make  possible 
investigations  of  the  quantitative  tradeoff 
between  optimal  performance  and  robustness, 
and  we  Illustrate  the  application  of  this 
differential  geometric  approach  via  various 
specific  examples. 

I .  Int  roduc  t ion 

It  is  well  known  that  there  is 
increasing  interest  in  the  employment  of 
robustness  techniques  for  the  discrete  time 
detection  of  signals  in  imperfectly  known 
noise.  The  traditional  approach  toward 
addressing  questions  within  this  rather 
broad  area  of  research  has  been  to  rely 
heavily  on  the  classical  saddlepoint 
criterion  of  Huber  (see,  for  example, [1] )  . 

A  variety  of  work  appearing  in  the  engi  - 
neering  literature  has  verified  that  such 
an  approach  can  lead  to  tractable  results. 
However,  it  may  be  argued  that  the  degree 
of  robustness  obtained  owes  much  to  the 
types  of  noise  models  admitted  by  the 
method.  In  reality  it  may  not  be  easy  to 
verify  that  the  types  of  models  appropriate 
to  the  saddlepoint  robustness  approach 
sufficiently  represent  the  full  extent  of 
variation  of  the  unknown  perturbation  of  a 
distribution  around  the  nominal.  Moreover, 


although  it  is  possible  via  the  approach  of 
[2]to  obtain  gener al r ep r e sen t a t ions  of  the 
noise  model  via  Choquet  capacities  [3],  it 
has  yet  to  be  seen  if  such  elegant  methods 
are  capable  of  enhancing  the  denseness  of 
the  class  of  noise  models  beyond  the 
relatively  few  standard  models  (see,  for 
example,  [1]).  In  addition,  the  saddle  - 
point  criterion  is  inherently  a  nonquanti- 
tative  approach  toward  imparting  robustness. 
We  intuitively  might  suspect  that  robust  - 
ness  is  obtained  by  a  judicious  tradeoff 
with  optimal  performance,  and  we  thus  might 
desire  a  way  to  quantitatively  measure  the 
degree  of  robustness  in  order  that  a 
weighted  combination  of  robustness  and  per¬ 
formance  could  be  considered  subject  to 
some  cost  criterion.  In  this  paper  we 
present  an  entirely  different  approach 
which  views  the  robustness  question  not 
from  the  saddlepoint  perspective  but  from 
one  which  is  rooted  in  differential 
geome  try. 


Development 


Viewing  the  robustness  problem  from  a 
slightly  different  perspective,  let  t>n 

denote  the  class  of  n-dimensional  distri¬ 
bution  functions.  From  this  point  of  view, 
the  performance  of  the  detector  is  thus 
expressed  by  considering  the  performance 
functional  we  then  simply  wish  to 

choose  the  detector  so  the  P  is  reasonably 
high  and  doesn't  vary  much  near  the  nominal 
element  of  £?n  .  Viewing  P  as  a  height 
function  over  ,  we  could  say  that  a 
robust  detector  would  yield  a  "surface" 
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which  above  the  nominal  element  is  both 
relatively  high  and  not  strongly  sloped. 

Such  a  perspective  thus  would  Indicate 
that  a  geometric  approach  to  robust  detect- 
tion  might  be  appropriate.  What  would  be 
needed  would  be  to  provide  a  differentiable 
structure  to  so  that  the  concept  of  slope 
would  have  the  proper  meaning.  We  would 
then  be  considering  a  height  function  over 
a  differentiable  manifold  M  which  would 
result  in  a  new  manifold  Mj  for  which  the 
Riemannian  metric  would  yield  a  norm. 

In  this  paper  we  present  some  specific 
applications  of  the  above  observations. 
Noting  that  a  Neyman-Pear son  approach  in  - 
volves  comparing  the  sample  vector  to  the 
appropriate  n-d imens ional  Borel  set  8n  in 
Rn  ,  where  n  is  the  number  of  samples,  we 
then  observe  that  in  this  robustness  appli¬ 
cation  we  could  in  practice  regard  8^  as 
specified  via  the  choice  of  nominal  distri¬ 
bution  under  Hq ;  we  then  would  be  interest¬ 
ed  in  analyzing  the  degree  of  variation  in 
the  false  alarm  probability  a  and/or  detec¬ 
tion  probability  8  as  the  under ly ing  d is t r i- 
bution  varies  about  the  nominal,  thus  fix  - 
ing  a  choice  of  height  function  h :  R  -►R  for 
some  natural  number  m,  where  h( ■ )corres  - 
ponds  to  the  value  of  a  or  8  for  some  fixed 
detector  of  interest. 

Consider  first  the  case  where  the 
class  of  n-d  imen s iona 1  distribution  func  - 
tions  is  parameterized  by  m  parameters; 
this  class  can  then  be  identified  with  a 
subset  of  R™ ,  and  the  corresponding  mani  - 
fold  Mjls  a  surface  in  Rm^.  An  appropriate 
metric  tensor  g(-,-)  is  inherited  from  the 
standard  Inner  product  on  Rm  ^  with  the 
obvious  choice  of  coordinate  system 


•y^  11,  12,...,  Ira,  'x^  leading  to  the 

components  of  the  metric  tensor  given  by 
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Associating  the  slope  of  the  unit  normal 


with  the  cosine  of  the  angle  of  the  unit 
normal  to  vertical  (with  Mj  immersed  in  R  +  ^ 
it  is  then  straightforward,  although  some¬ 
what  lengthy,  to  show  that  at  the  point 
corresponding  to  the  nominal  distribution 
this  cosine  is  given  by 


cos 


(1+  l  ( 

i-1 


2  -1/2 

)  ) 


Note  that  y  provides  a  measure  of  local 
m 

"first  order"  robustness;  smaller  values  of 

y  suggest  less  variation  in  a  or  6  near  the 
m 

nominal  distribution. 

Consider  now  the  discrete  time  detec¬ 
tion  of  a  constant  signal  s  in  additive 

i.i.d.  Gaussian  noise  with  mean  u  and 
2 

variance  o  .  We  note  that  there  may  be 

some  uncertaintly  in  all  of  the  values  of 
2 

s,u  and  a  .  Employing  first  the  linear 
detector,  we  then  choose  h(-)  to  correspond 
to  8  and  then  straightforwardly  obtain  (for 
n  samples) 

“  (n/(2’T02)),4exp(-[n  (  P  +  s  ) -T  ]  2  /  (  2nc?)  ) 


^•2  *•-  (  n  (  U  +  s)  -T)  • 

•exp(-[n(U+s)-T]2/ (2no2))  /  (Simo6)4  , 

where  the  threshold  T  is  specified  for  a 

given  false  alarm  rate  a  by  evaluating  the 

2 

detector  at  the  nominal  values  of  s ,  u,  ando 
We  next  employ  the  robustified  version  of 
the  linear  detector,  which  replaces  the 
identity  function  of  the  linear  detector 
with  the  nonlinearity  g(-)  defined  by 


g(x) 


if 

if 

if 


x  >  k2 

kll  x  i  k2 

x  <  k , 


It  is  well  known  that  this  detector  resists 
the  tractable  development  of  closed  form 
expressions  for  oiorS  in  the  Gaussian  case. 

In  order  to  numerically  compare  the  robust¬ 
ness  of  this  "censored"  version  of  the 
linear  detector  we  therefore  employ  a  large 
sample  Gaussian  approximation  of  the  test 
statistic.  The  resultant  lengthy  analysis 
shows,  for  example,  that  with  n*  50  ,  y*0 . 0  5, 
kj”-0.4,  k2”0.6  and  nominally  u*>0 . 1 ,  -  s*»0.  t ,( 

2  •  ■  o'  '  .  i 

'  -1,  we  have  y  ^  *34.. 1  ,.  .which  can  be  •  com-i 

pared  'fct)  -y  ^  *68 . 8°  for  the  linear  detector. 
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The  robustness  of  the  detector  which 
employs  censoring  is  thus  quantitatively 
demonstrated . 


Ill  . 


The  General  Case 


It  would  also  be  important  to  consider 
more  general  classes  of  distributions; 
ideally  we  would  wish  to  place  virtually 
no  constraints  on  the  admissable  distribu¬ 
tion  function.  Such  an  approach  is  actu  - 
ally  feasible  in  the  i.i.d.  case.  Since 
a  and  B  are  expressed  via  an  integral  over 


a  Borel  set  8  with  respect  to  the  appro- 
n 

priate  n-d imens iona 1  distribution  under 


Hq  and  respectively,  we  can  without  loss 


of  generality  note  the  independence  of  the 
observations  and  investigate  perturbations 
in  ci  and  B  by  limiting  consideration  to 
the  class  of  those  univariate  distribu  - 
tions  given  by  step  functions,  i.e.  those 
functions  of  form 


m+1 


F  (  •  )  -  Z  a,  I.  (')  ,  where  the  inter- 

i-0  1  Ai 


vals  A,^  partition  R  and  we  take  a^-O  and 


am+1”l.  For  a  fixed  finite  partition  P 


of  R,  we  note  that  the  corresponding  class 


of  step  functions  can  be  viewed  as  par 


ameterlzed  by  elements  of  R  .  Letting 
F(-)  denote  the  nominal  univariate  distri¬ 
bution  of  the  observations,  we  then  can 
employ  the  aforementioned  methods  to  obtain 


an  expression  for  cos  y  where  for  each 

m , 


partition  a  Stieltjes  approximation  to 
F(’)  is  chosen  and  regarded  as  nominal  for 
the  parameterized  case.  We  then  define 


•lim  .whenever  the  limit  exists.  The 
IH-O  m 


number  i  may  be  thus  be  interpreted  as  the 
angle  of  the  unit  normal  to  vertical  for 
the  general  ( non p a r ame t e r  i  z ed )  distribution 
case,  and  as  before,  we  would  like  this 
angle  to  be  small  for  robustness.  We  also 
have 


m  ~  2  -4 

cos >“lim  (1+  Z  ( — : )  ) 

|fl-0  i-1  13  i 


where,  as  before,  we  limit  consideration 
to  che  situation  where  the  height  function 
h(-)  is  i  or  (with  the  corresponding 
univariate  distribution  of  the  observa  - 


tions  F(-)).  Now  let 


3*  8n  “  { (x^ ,x^ ,... ,xn ): there  exists  £>0 


such  that 


(X1 . Xi-l*y*xi  +  l’-’-*xn)<  6n 


for  y  e  (x^-e,  x^)  and 


.i_l, z,Xi+i,...xn)eBn 


for  ze  (xi>xi+e)} 


3,  8  -  {  (x,  ,x.,...,x  ):there  exists  £>0 

in  1  c  n 


such  that 


(x1,...,x1_1,y,xi+1 . xn),Bn 


for  y  e  (x^E,  x  ) 


and 


(x1,....xi_1,z,x1+1,...,xn)t8n 


3H8n-  {  (xl . Xi-1 


for  z  t  (x  ,x  +e)  } 


l,xi+l . \i,yl . 


* .  .y. 


yj  _i  *yj  +  i . y n )  ;  there  exists  w  such  that 


(x  ,...,x  ,w,x ,X  )  f 3  8  , 

1  i-1  i  +  1  n  i  n  and 


(y  i»  •  •  •  »y  j-i*w ,y  j+i . yn}  e5j  8n} 


Similarly,  we  define  3^  by  replacing  3* 


in  the  3**  expression  with  3^.  In  an 


-+ 


analogous  manner,  we  also  define  3  ^  and 


3  .  We  then  can  establish  the  following 


result,  which  provides  a  closed  form  ex  - 
pression  for  cosy: 

Theorem:  Suppose  that  if 


xeR"«  3  |b  "3  ^B  *  ...*’3+8n3'B<'3,Bn  ...n3~8 
ln2n  nnln2n  nn 


then  x.«  int  (8  )uint  (?  )  .  We  then  have 
n  n 


2 

cosy*(l+A  +D  ,  where 


I, 


iBn 


dFCy^  .  ..dF(yi_l>dF(y1  +  1) 


•  d  F  (  y  n  )- 


i8n 


dF(Vl) . . ,dF(y1_1)dF(yi+1) . . .dF(yn) 


"-Z  / dF(x 

ii'j  + 

“  11 


) . . .dF(x1_1)dF(x1  +  1) . . *  d  ^ ( *  ) 


ij  n 


dF(y1)...dF(yj_1)dF(y1)...dF(yn)  + 


+  ydF(x1)...dF(x1_1)dF(x1  +  1)...dF(xn)dF(y1) 

' ij 8n  “•dF<yJ.i)dF(yj+1) • • -dF(yn) 


/ 


l' 


CiAur 

JNSPECTEO  j 


. . .dF(yj_1)dF(yj+1) . . .dFCy^)- 


-  /  dF(Xl)  .  .  .dF(x1_1)dF(xi  +  1)  .  . .dF(xn)dF(yi) 
ij  n 

...dF(y  . ) dF(y  , ) . . . dF(y  ) .whenever  the 

J “ 1  j+i  n 

integrals  exist. 

Proof :  Recall  the  step  function  approxi  - 


mation  F(-)  to  the  nominal  univariate 
_  m+1 

distribution  F(’),F(*)*  L  a  I  (•).  where 


i-0  1  Ai 


ag»0  and  am+^  “1.  The  associated  approxi¬ 
mate  measure  of  performance  is  given  by 

h(-)“Jg  dF  (y ^ )  .  .  . dF (y  ) ,  and  thus 
n 

m  m 


~V  :-r  nai(I8<xi-i-xi  ••••*!>- 

.?a,  i  *0  i  =0  n  2  n 

i  2  n 


~ls  (xi,xi  . xt  ))  + 


n  i  n 

mm  m  . 

+  -  £  I  a  _  ( I  (x  ,  x  ,  x  . x  ) 

Oij-  0  in-0  8n  h  1-1  i3  1n 


-Ig(Xi  .X  X i  .....X  )) 

n  1  3  n 


+  .  . ,+L . . . E  a  ( I  (x  ,  . 

i  -  0  i-0  n  8  h 

1  n-i  n 


1  i  .x.  J 

n-1  1-1 


. xi  ,xi^  )  * 

n  - 1 


8  '  i, 

n  1 


where  afc-(a  -a  )...(a  -a  ) 

1  1  k-1  k- 1 


<<,i.J.1+i“ai.)--,(ai+rai  >  and  xi’sup  Aifor 
k  +  1  k  +  1  n  n 

m 


1*0,1,.  ...m.  Expanding  the  E  ( — —  ex  - 


i  -1 


pression,  we  observe  that  as  the  norm  of 
the  partition  approaches  zero,  the  squared 
terms  collectively  reach  as  a  limit, 
whereas  the  cross  products  yield  T. 

QED 


We  remark  that  from  an  intuitive  per¬ 


spective,  consists  of  those  elements 


of  the  boundary  of  8^  which  are  intersected 


by  rays  parallel  to  the  i  axis  and  moving 
in  the  positive  direction  from  the  exterior 


of  8  to  its  interior,  whereas 
n 


3. 8  is 
1  n 


-  («<*!>  •  •  •^(xi-1)dF(x1+1)  .  .  .dF(xn)dF(ji) 
"1J  n 


formed  in  an  analogous  manner  with  the  rays 
moving  in  the  negative  direction  from  the 


exterior  of  8  to  its  interior.  In  addition, 
n 


we  note  that  the  Theorem's  hypothesis  is 
simply  a  mild  condition  on  the  regularity 


of  8^  which  is  frequently  easy  to  satisfy 


in  this  detection  context  (wherein  8  arises 

n 


by  way  of  a  threshold  comparator).  Moreover, 
the  existence  of  the  integrals  in  the 
Theorem  is  very  often  easy  to  verify  since 


the  boundary  of  8^  is  sufficiently  well  be¬ 


haved  in  such  cases. 

As  an  example  of  an  application  of  the 
Theorem,  consider  again  the  1 inea r de t e c t o r . 
For  n-2  it  then  follows  that  (where  the 
detector  threshold  is  T) 

Sl82”  32B2-  Hx.y)  :y“T-x}  ,B~82»  B 2 - 


.++R  , 

^l  2  82 


^2 T 8 2 ”  ^x>y)  ;y~x } 


31282~  32  1  B 2 “ 


31282* 


32182* 


^1282“  h  182‘ 


<t>. 


We  therfore  have 

cosY-l/(l  +  l  +  l+0)!i-l/3,s,  i.e.  y-54.70, 
gardless  of  the  nominal  distribution, 
may  be  generalized  to  show  that  for  n 
pies 

_L 

cosy*(l+n)  ,  regardless  of  the  nominal 
d is t r ibut ion  . 


r  e- 
This 
sam- 


Note  that  lim  cosy-0, i.e. 

n-*-oo 

the  linear  detector  becomes  completed  un- 
robust  (as  measured  by  y  )  as  the  number  of 
samples  approaches  infinity. 

On  the  other  hand,  consider  the  class¬ 
ical  robustifled  version  of  the  linear  de  - 
tector,  wherein  a "censored"  detector  non¬ 
linearity  g(')  of  form 


[  x 

if 

1*1 

<k 

gU)- 

k 

if 

X 

>k 

is  used 

1  -k 

if 

X 

<  -k 

It  then  can  be  shown  that 


1 B^m  { (x , y ) : x-T-k, y >k ;  or  y«T-x, 


T-k<x  <k } 


,B,“{  (x,y)  :y-T-k,  x'-k;  or  y-T-x, 

T-k  <x  <k } 


1182- 


^2  82 " 


3++8  • 
1  2  °2 


+  +r 


>2  ’  83“  {  ( x  ,  y  )  :  x  and  y  >k  ,  or  y-x  and 


T-k  <x  <k 


R+  ■  8  -  > 

12  2  21  2 


*  1 2  82* 


r+8  ■ 

212 


1  2 82*  ’2 1 82* 


For  this  case  the  exact  value  of  cosy 


depends  on  the  values  of  k,T,  and  the 

choice  of  nominal  distribution.  However, 

we  can  make  some  general  conclusions  when 

the  amount  of  censoring  a pproache s  maxima  1 

(k-0)  .  In  this  case  we  have,  when  n«2  , 

lim  cosy  -(1  +  2(1-F(0))2  +  2(1-F(0))2)_is. 
k-0 

For  the  common  case  where  F(0)^  we  then 
obtain 

lim  cosy  <_  1/2^  ,i.e.  limy^  45  ,  which 
k-0  k-0 

may  be  compared  to  the  linear  detector's 
r»54.7°  for  n«2.  This  may  be  generalized 
through  a  lengthy  analysis  to  conclude 
that  for  n  samples, 

lim  cosy  » ( 1+n ( 1 - F ( 0 ) ) 2 ^ n  ^+n(n-l)* 
k-0  , 

•  (1-F(0))n)‘\ 

Note  that  for  F(0)>0  we  have  lim(lim  cosy) 

n  — “  k  -*-0 

•  1,  that  is,  the  detector  approaches 

possessing  complete  robustness  (as  measured 
by  ,  )  as  the  number  of  samples  tends  to 
infinity.  For  F(0)<^  ^  we  have 

lim  cosy  <(l+n-2"2(n-1)+n(n-l)2'n)"S5. 
k-0 

Note  that  the  upper  bound  can  be  approach¬ 
ed  arbitrarily  closely  for  F(0)  near  4 • 

For  n-3  this  upper  bound  becomes  0.8 
(corresponding  to  >”36.9°),  whereas  for 
n»10  it  becomes  0.96  (corresponding  to  y 
”16.3°) .  This  can  be  compared  to  the  case 
of  the  linear  detector,  where  for  n"3  we 
have  »60°  and  for  n»10  we  have  y“72.5°. 

For  the  larger  values  of  n  the  robustness 
advantages  of  the  classical  robustified 
linear  detector  are  thus  quantitatively 
demonst  rated . 

IV.  Second  Order  Robustness 

Finally,  we  note  that  our  geometric 
approach  admits  the  additional  "second 
order"  sensitivity  check  provided  by  cur¬ 
vature.  There  are  many  different  defini  - 
tions  of  curvature,  however  scalar  curva¬ 
ture  has  the  advantage  of  being  indepen¬ 
dent  of  the  local  coordinate  system  em  - 
ployed,  thus  simplifying  its  computation. 

In  addition,  the  choiceof  scalar  curva¬ 
ture  is  intuitively  appealing  since  it  is 
just  a  sum  of  the  various  sectional  curva¬ 


tures;  we  are  therefore  simply  accumulating 
ordinary  two  dimensional  Gaussian  curvature 
it  all  possible  orthogonal  directions.  For 
the  parameterized  distribution  case  it 
would  then  be  possible  to  generate  numer  - 
ical  values  for  the  "second  order"  measure 
of  robustness  provided  by  scalar  curvature, 
which  from  [4]  is  given  by  (where  here  the 

Einstein  summation  convention  is  used) 
l  0 

ikr3rik 

V8  [irr 


sr 


_ Li  +  r2-  ps 

3x,  s2  ik 


-  r 


sk 


rs  i 

i  l 


(m  not  summed),  where  the  Christoff  el 
k 

symbols  T  are  given  by 
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which  (g1^)  “  (g 


ij 


3x  . 
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)-1.  Although  the  above 


-  -J 

equations  regarding  scalar  curvature  appear 

rather  compact,  numerical  ca  1  cu  1  a  t  ioiis  in¬ 
volving  them  can  be  quite  tedious  (but  not 
difficult),  especially  for  large  m. 

III.  Conclusion 

We  have  presented  a  new  approach  to¬ 
ward  robust  signal  detection  which  is 
based  on  differential  geometric  methods  as 
opposed  to  classical  saddlepoint  criteria. 
These  techniques  are  seen  to  admit  a  quan¬ 
titative  measure  of  robustness  through  the 
geometric  concepts  of  unit  normal  slope 
and  scalar  curvature,  thus  allowing  the 
consideration  of  a  weighted  combination  of 
performance,  first  order  robustness  (via 
unit  normal  slope),  and  second  order  ro¬ 
bustness  (via  scalar  curvature)  subject  to 
some  cost  criterion  of  interest.  Our 
techniques  are  additionally  illustrated  in 
the  paper  through  various  specific  exam¬ 
ples. 
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